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Abstract 

This is paper I of a series of two papers, offering a self-contained anal- 
ysis of the role of electromagnetic stress-energy-momentum tensors in the 
classical description of continuous polarizable perfectly insulating media. 
While acknowledging the primary role played by the total stress-energy- 
momentum tensor on spacetime we argue that it is meaningful and useful 
in the context of covariant constitutive theory to assign preferred status 
to particular parts of this total tensor, when defined with respect to a 
particular splitting. The relevance of tensors, associated with the electro- 
magnetic fields that appear in Maxwell's equations for polarizable media, 
to the forces and torques that they induce has been a matter of some 
debate since Minkowski, Einstein & Laub and Abraham considered these 
issues over a century ago. The notion of a force density that arises from 
the divergence of these tensors is strictly defined relative to some inertial 
property of the medium. Consistency with the laws of Newtonian con- 
tinuum mechanics demands that the total force density on any element 
of a medium be proportional to the local linear acceleration field of that 
element in an inertial frame and must also arise as part of the divergence 
of the total stress-energy-momentum tensor. The fact that, unlike the 
tensor proposed by Minkowski, the divergence of the Abraham tensor de- 
pends explicitly on the local acceleration field of the medium as well as 
the electromagnetic field, sets it apart from many other terms in the total 
stress-energy-momentum tensor for a medium. 

In this paper we explore how electromagnetic forces or torques on mov- 
ing media can be defined covariantly in terms of a particular 3-form on 
those spacetimes that exhibit particular Killing symmetries. It is shown 
how the drive-forms associated with translational Killing vector fields lead 
to explicit expressions for the electromagnetic force densities in stationary 
media subject to the Minkowski constitutive relations and these are com- 
pared with other models involving polarizable media in electromagnetic 
fields that have been considered in the recent literature. 

1 Introduction 

The interaction of matter with the electromagnetic field has played a dominant 
role in the development of our understanding of Nature. In classical Newtonian 
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continuum mechanics one is concerned with dynamic processes involving the 
interaction of (bounded) classical continua with external forces and torques in 
space (see e.g. [T]). In relativistic continuum mechanics the theory is generalized 
to incorporate the concept of the spacetime manifold and formulated in terms 
of tensors on this manifold (see e.g. |2]). In both formulations the histories of 
configuration variables must be compatible with balance laws associated with 
possible translational and rotational symmetries of a metric structure and an 
associated balance law of 'energy'. If electromagnetic fields are involved these 
laws are supplemented with the macroscopic Maxwell equations in media. As 
a result of cohesive forces originating at the molecular level the interaction of 
a material with these (and gravitational influences) results in locally induced 
strains. These strains determine material stresses that are encoded into various 
'stress-tensors'. These in turn describe the local distribution of force and torque 
densities in the medium. The mathematical structure of such tensors charac- 
terizes the response of different materials to external influences. Constitutive 
relations are relations that are used together with field equations, equations of 
motion and boundary conditions to fix the dynamic evolution of the indepen- 
dent dynamical variables. If the configuration of the system (open or closed) 
involves thermodynamic variables the classical laws of thermodynamics can be 
used to constrain these relations for real media. In the Newtonian formula- 
tion one uses the Maxwell and Cauchy stress tensors [3], together with certain 
equations of state. In the relativistic formulation a primary role is played by 
the total stress-energy-momentum tensor for the system. The determination of 
these tensors, together with appropriate constitutive relations for different types 
of polarizable uncharged media often requires input from experiment. The elec- 
tromagnetic constitutive properties and the appropriate electromagnetic stress- 
energy-momentum for light in moving media has been a subject of debate (and 
possible confusion) for over a century. Difficulties arise in properly accounting 
for the local nature of the wave-matter interaction in terms of experimentally 
measurable effects. Even for static electromagnetic fields there are (often non- 
linear) subtle interactions that induce changes in shape or volume in deformable 
media dependent on its thermodynamic state. 

There is a considerable body of opinion suggesting that in some sense the 
choice between different stress-energy-momentum tensors describing interac- 
tions of a material medium with an electromagnetic field is a matter of con- 
venience and that different choices simply provide alternative descriptions of 
the same overall interactioi|3 [H |5l |6l I3 • Since there is no preferred tensor par- 
tition of the total stress-energy-momentum tensor into sub-tensors describing 
the behaviour of interacting sub-systems there can be no unique definition of 
a stress-energy-momentum tensor describing an interacting subsystem. Thus 
for electromagnetic fields interacting with an electrically neutral bounded con- 
tinuum it is always possible to redefine integrated electromagnetic forces and 
torques on the medium associated with sub-tensors according to taste, partic- 
ularly if the medium is composed of piecewise inhomogeneous material subsys- 
tems or the fields are time dependent. However in order to model the total 
interacting system the choice of a total stress-energy-momentum tensor is nec- 
essary. If one has decided how to model a neutral but polarizable medium in 

^ The genesis of this idea of classical duality may have its origins in the analogy with 
wave-particle duality in quantum mechanics. 
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the absence of externally applied electromagnetic fields then different choices 
of a stress-energy-momentum tensor for electromagnetic fields in the medium 
to describe the additional interactions with such fields will inevitable lead in 
general to different predictions for the interacting system. 

Although the laws of classical electromagnetism in the vacuum are firmly 
established in a relativistic context there is no general agreement on how best 
to accommodate the dynamics of deformable media as a self consistent theory 
on spacetime, (see e.g. [5]). This makes any rigorous formulation of relativistic 
continuum mechanics of inhomogeneous dispersive polarizable media interacting 
with electromagnetic and gravitational fields difficult even if one contemplates 
using it for systems in non-relativistic motion in some frame of reference. 

Despite these shortcomings there has recently been a resurgence of interest 
in the so-called Abraham-Minkowski controversy and its relevance to the use of 
either the Abraham or the Minkowski form of stress-energy- momentum tensor in 
interpreting experiments involving electromagnetic fields in media [Ul ll0l[Tnil2| . 
Part of this difficulty is no doubt due to the complex nature of material responses 
to forces in general. From our perspective such experiments are seeking con- 
stitutive relations involving particular material systems and the macroscopic 
Maxwell fields in media. As such it should come as no surprise that different 
systems might yield different responses particularly if the competing effects of 
electro- or magneto-striction mentioned above are contributing differently in 
different experiments. 

A number of historic experiments have sought to detect the discriminating 
Abraham force [131 114 | . This is strictly discriminatory only for homogeneous 
non-dispersive stationary media, so calls into question those experiments that 
involve media in motion. Indeed the fundamental difference between the forces 
or torques induced by the divergence of the Abraham and Minkowski electro- 
magnetic stress-energy-momentum tensors is that the former, unlike the latter, 
can depend explicitly on the acceleration of the medium. From this observa- 
tion it is our contention that, for any specified electromagnetic constitutive 
relation, there remain experimental avenues offering new means to discriminate 
between alternative proposals for the form of the electromagnetic stress-energy- 
momentum tensor in media, despite the attendant inherent material constitutive 
complexities. In particular we argue that, in principle, the observed dependence 
of a time-averaged electromagnetic wave-induced torque on the angular speed 
of a materially isotropic but inhomogeneous uniformly rotating electrically neu- 
tral medium could be used to discriminate between the electromagnetic wave 
interactions described by the Abraham tensor from those described by the sym- 
metrized version of the tensor introduced by Minkowski, and possibly those 
proposed by others. 

It is therefore of interest to calculate how such a torque depends on geometric 
properties of a cylindrical insulator, its speed of rotation and its electromagnetic 
constitutive properties. Since the medium will be assumed electrically polariz- 
able and magnetizable (but non-conducting) one is immediately confronted with 
a number of subtleties associated with the form of this tensor and the material 
constitutive properties of the medium. These issues have a bearing on how one 
formulates the classical electromagnetic force on a macroscopic body particu- 
larly one that is accelerating. Since this has led to a number of related questions 
in the recent literature (THl HH [H] , this article attempts to make explicit our 
perspective. After reviewing approaches in the Newtonian framework the use of 
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a fully covariant special relativistic framework is advocated. The essentials of 
this formulation are described below which should be read in conjunction with 
the expository material on electromagnetic theory in the language of differen- 
tial forms contained in the appendix. In the following paper, applications to 
rotating media are discussed in this context. 



2 Points of Departure 

In order to motivate our method of analysis leading to a computation of the 
electromagnetic torque on an inhomogeneous rotating uncharged insulator, it is 
useful to place our methodology in the context of the recent literature in this 
subject. Some of this literature is devoted to the derivation of expressions for 
the classical force (and torque) induced by electromagnetic fields on electrically 
neutral polarizable continua based on an underlying discrete model of electro- 
magnetic sources. A traditional non-relativistic approach is to take as a point 
of departure the classical vacuum Maxwell equations for the fields e^{r,t) and 
h^{r,t) in some inertial frame (here labelled U) and moving point source^ to- 
gether with the Newtonian equations of motion in for the sources. For a 
collection of N point charges where the a-th point has mass and Newtonian 
velocity Vait) at time t, the motion of each particle is given by a solution to the 
N coupled ordinary differential equations: 

d 



dt 



(m„v„(i)) = f^{t) + J2j'c.pit), a,P = l...N (1) 



where J-ap{t) is the inter-particle force and J- ait) the resultant force on the 
a-th particle due to all other influences. The non-relativistic interaction of 
the ambient Maxwell fields e*^ and b'^ with any point particle with charge Qa 
located at r = (t) with Newtonian velocity Vq, (t) = Tq, (t) is derived from the 
electromagnetic force -^^*^(i) on that particle with: 

= gc(e^(ra(t),t)+v„(t)xb^(r„(i),t)) (2) 

in the Gibbs vector notation. The fields themselves are in turn derived from 
the vacuum microscopic Maxwell equation with singular sources. In the Gibbs 
vector field notation these are the equations 

curle^ = (3) 
divb^ = (4) 

dt 

endive^ = (6) 



/iQ-'curlb^ = jy + e„^ (5) 



respectively, where A*o ^ — "^o^o, cq denotes the speed of light in vacuo, 'fP = 
J2a <laS^{r-ra{t)) and = J2a1a^ait)S'^{r-ra{t)) in terms of the singular 
Dirac distribution 6^ . In these equations only the fundamental electromagnetic 



We assume throughout that free magnetic charge is absent in Nature. 
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fields e*^ and play a role. Multi-pole magnetization and polarization sources 
are defined in terms of the attributes of the charged particles and their motion 
as limits in a multi-pole expansion. In non-relativistic models retardation effects 
due to the Maxwell displacement current are often ignored. It is straightforward 
to verify that the total linear particle momentum in the system, rnayait), 
is a constant of the motion provided J2a -^ait) = and J-ap{t) = —J-fja{t)- 
From these fundamental assumptions one may approach a continuum descrip- 
tion from a number of distinct directions including a multi-pole expansion of 
the electromagnetic fields about some arbitrary point in space followed by a 
spatial smoothing procedure for the singular sources, in order to generate a bal- 
ance law for smoothed out total non-relativistic linear momentum for the total 
interacting system. This yields an Euler continuum description from the funda- 
mental particle-field description. While the resulting overall continuum balance 
equation may not be sensitive to the precise nature of the smoothing functions 
the interpretation of individual forces in the Eulerian balance relation may be 

HOI [21]. 

An alternative approach to a continuum description has been to assume 
that the point particles are electrically neutral (atoms or molecules) but are 
endowed with elementary electric dipole moments and/or magnetic dipoles or 
elementary current loops. This requires that ^ be changed to refiect this 
modification. Subsequent smoothing would lead to a balance law involving force 
terms different from the point charge model even in the static limit. This should 
come as no surprise since the underlying microscopic models are different. 

In comparing the predictions of different continuum models one must be 
aware that they may differ only in their effects at the boundaries of spatially 
compact media. Terms that may be discarded during an integration by parts 
in the development of the modelling process may well contribute to boundary 
interactions depending on the nature of the boundary conditions. 

To illustrate some of these points consider the typical modelling of a homo- 
geneous dispersion-free electrically neutral stationary medium [22 based on a 
collection of point charges yielding, in some continuum limit, a time dependent 
volume force density on given by 

Ft = iv(e^.p^ + b^.m^) + i|(p^xb^) = F,,i + F.d + F,,k , (7) 

expressed in terms of time dependent electromagnetic vector fields on in a 
global orthonormal Cartesian frame with basis fc} and Cartesian coordi- 
nates {x^ , , x^} . The force vector on a finite volume Vt of the medium is 
then 

Ft = i Ft dx^ dx'^ dx^ 
Jvt 

In ([7]), the p^ and are spatially smoothed time-dependent electric and 
magnetic dipole fields on R"^ obtained by truncating a particle force multi-pole 
expansion. Prior to smoothing they can be expressed in terms of the charge, 
position and velocity of the electrically charged constituents of the medium 
relative to some arbitrary point and the bulk velocity in the medium. They are 
superscripted to indicate that these fields are referred to an inertial (laboratory) 
reference frame U . In principle p'^ and raP are determined from the continuum 
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limit of the particle equations of motion. In practice this is difficult so one 
resorts to constitutive relations and assumes a form for the bulk motion of the 
continuum. With these closure relations the vacuum Maxwell system absorbs 
the smoothed polarization and magnetization sources into the phenomenological 
fields p*^, m'^, d*^ and h'^. 

In view of the covariant language used in the rest of this paper, it is useful to 
write ([7]) in tensor form using the Killing symmetry of the Euclidean structure 
defined by the metric of 3-dimensional space. The reader should consult the 
appendixfor a self-contained formulation of the Maxwell system in terms of 
the frame dependent differential 1-forms e*^, b*^, d^, h^, p*^, on space and 
their definition in terms of a unit time-like vector field U and the Maxwell and 
excitation 2-forms F and G respectively on 4-dimensional spacetime endowed 
with the Lorentzian metric tensor field (jj. In a Cartesian coordinate system on 
M'^, the covariant and contra- variant metric tensor fields are 

g - ^da;*^da;\ with g-^ = J^A-^l-. 

i—l i—1 

and in these coordinates the associated translational (Killing) vector fields are 
{jfi:} (for all time t) satisfying £_g_g = 0. If is any of these Killing fields, 

theifl K = dx*, for some Cartesian coordinate x*, i = 1,2,3. In terms of 
K = with #1 = dx^ A dx^ A dx^, a Euclidean covariant representation of 
([7]) that meshes with the conventions to be established below is given in terms 
of the equivalent force 3-form for direction K: 

:Fk = i[d(e^-p^ + b^.m^)+#£a.(p^Ab^)] A^, (8) 

where we define a ■ /S = (a, /3) = g ""^ (a, /3) for any spatial 1-forms a, /3 and 
d denotes the spatial exterior derivativq^ on M'^. Then the Cartesian compo- 
nent of the vector force in the Cartesian direction i on a volume Vt due to its 
surroundings and ambient fields in the medium is given by 

JVt 

Note that in the static field situation the force 3-form ([75]) is expressed in terms 
of the gradient of the interaction scalar ^(e'^- p'^ -I- b'-^- m'^). Using the macro- 

^ All electromagnetic tensors in this article have dimensions constructed from the SI di- 
mensions [M], [L], [T], [Q] where [Q] has the unit of the Coulomb in this system. We adopt 
[g] = [L^] , [G] = [j] = IQ] , [F] = 1^ where the permittivity of free space eo has the dimen- 
sions [Q'^T^ L^^]. Note that the operators d and V defined in the appendixpreserve the 
physical dimensions of tensor fields but with [g] = [L^], for p-forms a in 4 dimensions, one 
has [*a] = [a][L4-2p]. 

*For any metric tensor field Q, we define the 5-dual of any 1-form a to be the vector field 
5 = Q~^(a, —). Similarly, for any vector V, we define the 6-dual 1-form V = Q{V, —). 

^ See appendixfor the definition of the spatial exterior derivative d in terms of the exterior 
derivative d on spacetime, and for definitions of other spatial operators. 
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scopic Maxwell equations 

de^ = (9) 

dB^ = (10) 

dh" = + -b^ (11) 

dB" = (12) 

where D'-^ = #d'^, B*^ = #b^, a = coCg^a for any p-form a and with the 
definitions 

J" 1 / • \ 

-^ = -ijjjp = — (pU^dm'^\ and p';! = -{iu * jp) * U = -dP^ , (13) 
Co Co V / 

one may write0, with = # =F and E'^ = #e^: 



p^e^{K) + i^jl^ A K + C^m.^ A K + CkP^ A-E^ + d{e^ iK)P^) . (14) 



Although the overall charge of the medium is taken to be zero the first two terms 
in the equation for J-fc above constitute a local Lorentz force densitjH generated 
by the induced polarization charge p^ and current ■ The exact 3-form in the 
last term will contribute to forces on the spatial boundary of the medium unless 
they happen to be zero as a result of boundary conditions satisfied by the fields 
there. The remaining terms describe local forces depending on inhomogeneities 
arising from the spatial rate of change of polarization and magnetization in the 
system. 

The above derivation of the structure of a local Newtonian force density 
in a neutral macroscopic continuum in an electromagnetic field starts with a 
particular non-relativistic model in an inertial frame. We stress that in the ab- 
sence of dynamical information about the polarization and magnetization the 
cogency of (1141) demands a knowledge of supplementary constitutive relations to 
accommodate the response of the medium (via p'^ and m'^) to the fundamental 
fields e'^ and b'^. One could proceed to express this force density as the sum 
of the spatial divergence of a time dependent Maxwell second rank stress ten- 
sor on and a time derivative modulo boundary termtj^]. Adding the Cauchy 
tensor describing the medium in the absence of electromagnetic forces to such 
a Maxwell stress tensor would give the total Cauchy stress for the medium (in 
the presence of fields) that enters into the local balance law for non-relativistic 
linear momentum for the complete interacting system. The classical local dy- 
namics of this system follows as a solution to such a balance law and requires 
implementation of interface conditions at media boundaries or interfaces where 
material properties of the continuum are discontinuous. The latter constraints 
benefit from a distributional reformulation that has been discussed more fully 
elsewhere [H 123111]. 

An alternative point of departure for the modelling process is to start with 
a fully covariant total symmetric stress-energy-momentum tensor T*°* for a 
medium interacting with classical fields on a general spacetim^l. The fully 



^See appendix for the definition of the Hodge maps # and *, and further details of this 
calculation. 

^The factor ^ arises since the polarization is smoothly distributed in the medium. 
^ There is no unique way to perform such a split. 

®In the SI system the tensor T'°* has the physical dimension of a force, i.e. "'^Z' . 
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covariant local classical equations of motion of the medium are then postulated 
to be given by 

v-r*°* = (15) 

together with equations for the fields. Here the divergence operator V- is defined 
with respect to the Levi-Civita connection V on spacetime. The fundamental 
property of T*°* is that it acts as a source of Einstein gravitation. In most clas- 
sical considerations one demands T*°*([/, t7) > for all future pointing timelike 
unit vector fields U, reflecting the attractive nature of gravitation and T*°*(U, U) 
is then identified with local mass-energy density in the frame U . In a source- 
free region of spacetime containing only electromagnetic fields this condition is 
maintained and identifies field energy. 

In a spacetime with local isometries generated by a set of Killing vector 
fields associated with the spacetime metric g, r*°* can be used to generate a 
set of closed 3-forms on spacetime. The vanishing of the exterior derivative of 
each form in this set can in turn give rise to a conservation law when integrated 
over a regular 4-dimensional domain of spacetime provided the forms are free of 
singularities there. Then a class of Killing vectors that generate spatial trans- 
lations can be used to construct conservation laws that reduce to the balance 
laws for the components of Newtonian linear momentum in some non-relativistic 
Newtonian limit. 

Thus if T is any symmetric rank 2 symmetric tensor on spacetime so Tab — 
Tba where 

in any cobasis of 1-forms {e"^} with a dual basis of vector fields {Xi,} one defines 
the drive 3-form associated with T and K: 

= -CKT{K,Xa)*e'' (16) 

for any vector field K on spacetime where C,k = ±1 will be defined below 
to conform with the physical interpretation of the different components of the 
drive-form. Given a frame defined by the unit timelike (future pointing) observer 
fiel41 U, one may decompose into a spatial 2-form and spatial 3-form 
relative to U on spacetime: 

Tk = J^AU + p'^ (17) 

with i[/J^ — iupSi — 0- When K is a, Killing vector field {CkQ ~ 0) it is a 
mathematical identity (Benn & Tucker 1988; Benn 1982) thalQ 

dT^ = -CA(V-T)(if)*l. (18) 

Hence, if V • T = then 

dTji = 0. (19) 

"The frame is inertial if VU = 0. 

'^More generally if T is an arbitrary (2, 0) tensor, T = T^b e°'®e^ , of no particular symmetry, 
one has for any vector field W the identity ^{Cwg){^a,Xf,)T^°-^^ ★ 1 = — Cvk^^tiv — (V • 
(SymT))(W) * 1 where —C,w ''W = {SymT){W, Xa) * and SymT is the symmetric part of 
T with components T^^^y = i(T„i, -|- T^a)- 
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In terms of Jj^ the conservation equation (|19p becomes 

dJ^ + CuTj^ = 0. 

If if is a spacelike translational Killing vector field with open integral curves 
then 

is a linear momentum current (stress) 2-form in the frame U and 

is the associated linear momentum density 3-form in the frame U . If K is 
a spacelike rotational Killing vector field generating 50(3) group isometrics 
with closed integral curves, then jj^ is an angular-momentum current (torque 
stress) 2-form and is the associated angular-momentum density 3-form in 
the frame U . If X is a timelike translational Killing vector field, then is 
an energy current (power) 2-form and is the associated energy density 3- 
form in the inertial frame field U . In the following, attention will be directed 
mainly to particular translational and rotational spacelike Killing vectors K of 
fiat spacetime and the computation of integrals of j]^ and for a particular 
contribution to associated with electromagnetic fields in media undergoing 
various states of motion observed in an inertial frame defined by U . 

As noted above when T*°* describes a domain of vacuum spacetime, free of 
matter but containing electromagnetic fields, one requires the electromagnetic 
field energy in any local frame U to be positive, i.e. > in terms of 

p^j. Furthermore we shall require that, for any spacelike Killing vector field K , 
4fp^ > 0. This ensures that when K generates spatial translations in Minkowski 
spacetime and U defines an inertial frame then the time-averaged integral of pl- 
over a finite region of space can be identified with the time-averaged component 
of physical linear momentum associated with a harmonic plane wave in the 
direction of its propagation. These conditions are ensured if 

_ g{K,K) 
15(^,^)1 

i.e. C,K = 1 if if is spacelike and C/f = — 1 if if is timelike. 

If the spacetime admits a foliation by hypersurfaces t = constant then (fT9|) 
takes the form adapted to the frame U: 

dJ^-V^ = 0. 
Co 

In this spacetime framework, a basic postulate is that the history of a ma- 
terial continuum interacting with gravity and electromagnetic fields can be de- 
termined from and the Maxwell system (HTl) 

dF ^ Q and = j, (21) 

where the excitation 2-form G depends on the interaction with the medium and 
the 3-form electric 4-current j encodes the electric charge and current source. 
Such an electric 4-current describes both (mobile) electric charge and effective 
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(Ohmic) currents in a conducting medium. To close this system in a background 
gravitational field, electromagnetic constitutive relations relating G and j to F 
are necessary. These relations may also depend upon properties of the medium, 
including its state of motion. 

If one makes the arbitrary split T*°* = j-matter^BAf ^ j^em ^^^^^ general 
y ,j.matter,EM ^ _y.'j^EM ^ q ^ nou-relativistic modeUing situation, 

any interpretation of various terms in the decomposition of V • T*°* must be 
understood to be with respect to a particular splitting. This is of particular 
relevance in situations where V • T^^^ contains a coupling of electromagnetic 
fields and mass-energy to the bulk local acceleration field of a medium. 

In a fully coupled system that consistently incorporates the gravitational 
interactions using Einstein's gravitational field equations with a symmetric Ein- 
stein tensor, ((15]) becomes an identity. In descriptions with non-dynamic gravi- 
tation (jl5[) is part of the coupled system for the remaining dynamic fields. Since 
the effects of gravity will be neglected in the following we assume henceforth a 
background Minkowski spacetime and all observers will be inertial with VC/ — 0. 

Consider a material continuum macroscopic model in which r*°* contains 
the symmetric electromagnetic stress-energy-momentum tensor 

T^*^ = ^(^iaF® i^-G + iaG (g) ^F + ★(F A ★G).?) . (22) 

The excitation 2-form G in this expression must be specified in terms of F and 
other properties of the medium. For a simple non-dispersive isotropic medium, 
one has the constitutive relation 

G = eoerivF AV - — ir (iy * F Av) ^ ( er - —] iyF AV+ —F, (23) 

fir ^ / \ fJ-r J fJ-r 

where ^ is a unit, timelike 4- velocity field describing the bulk motion of the 
medium and and fir are relative permittivity and permeability scalars on 
spacetime. For a spatially inhomogeneous medium these will not be constants: 
der 0, d fir ^ 0. Furthermore, if the medium is accelerating then V\/F ^ 0. 
For any Killing vector field K, the tensor gives rise to the Killing drive 

3-form: 

t|*' - -^(iKGA^F-FAiK^G), (24) 

where the forms F and G are required to satisfy the Maxwell equations (H71) . 
For an uncharged non-conducting medium j — 0. The stress-energy-momentum 
tensor l|22p is that obtained by symmetrizing the one proposed by Minkowski 
to describe electromagnetic stresses and energy balance in a medium. 

3 Covariant Forces and Torques in Media 

In general the total stress-energy-momcntum for a medium may contain singu- 
larities and discontinuities in its material properties. To facilitate the following 
it proves convenient to restrict to a bounded domain of spacetime containing 
(the smooth history of) the medium immersed in the vacuum (see figure [ij. 
Thus we explicitly leave out of the discussion sources of stress that may arise 
from singularities or discontinuities in the medium history. One may regard the 
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Figure 1: Immersion of a material history in spacetime 



immersion as the material body and it is assumed here to have a topologically 
trivial structure with a total interior Killing drive form r]^ with interior sup- 
port including the smooth boundary of the body. The exterior of the body in 
spacetime is assigned a total Killingfl drive t]^"* with exterior support excluding 
the boundary of the body. 

For any observer field U, spacetime domain M and total A'-drive on M., 
one has from (fT7)) 

Jk = -t-uTk and = -[lu^rj^ )*U. 

If one writes 

= -CKclp„,e9iV,K)*V + T^ 

to describe a medium with bulk 4-velocity V on Ai and proper inertial mass- 
energy density scalaid Pme, given the conservation law d{pme * ^) = 0, the 
equation of motion dr-^ = becomes 

C,Kclprr.eA{K) - /j^ , 

where = — * dr^^ and the 4-acceleration vector field of is A = VyV . 
Contracting the * Hodge dual of this equation of motion with U and integrat- 
ing the resulting 3-form over any U-orthogonal spacelike hypersurface Et that 

*When the effects of gravity are neglected the interior and exterior spacetimes admit the 
same set of KiUing vector fields. 

^The first term on the right in the above equation arises from a contribution Cg pmeV Cg) V 
to the total stress-energy-momentum tensor. 
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intersects the domain M yields 

where the mass-energy 3-form = —CQPme * U and the total instantaneous 
integrated X-drive on Sf at time t in the U frame is 

Thus for any part t^'-' of the total iC-drive — J^j '^li''' ^ call the 
3-form —ijj dr'^''' the j-th part of the instantaneous X-drive density 3-form. 

Using this notation consider the two regions Mi = in and M2 = out with 
t'^ on A^i describing the interaction of a medium with electromagnetic fields 
according to d r]^ = and rj^"' describing electromagnetic fields in the vacuum 
according to drjp* = on A^2- The interface between these regions will be 
denoted S12 regarded as an immersion in spacetime. In the following the dif- 
ference of J^'™ and J^°^* across S12 will be encountered. This discontinuity 
jump is defined in terms of the pull-back map these forms to E12: 

A/ jU \ v'* f jU,in jU,out\ 

12[Jk) — ^12\-'k ^ -^K )■ 

With a split given by j = {matter, EM} write 

in in, matter , in, EM 

where j-^-'"'^**'^'" describes matter without permanent polarization or magneti- 
zation in the absence of external electromagnetic fields, the equation drJJ' = 
yields 

rU, in, matter i . rU,in.EM r-^ i /'oir\ 

Ik [^t\ + Ik ' l^t] = 0. (25) 

As noted in the introduction the electromagnetic constitutive relation may imply 
a coupling of the electromagnetic fields in the medium to its deformation tensor 
and hence may contribute to its strains. For an uncharged medium, if there are 
no electromagnetic fields to polarize the medium, the second term on the left 
in (j25p will be zero. In this case one may identify it with the A'-drive on the 
medium due to an applied electromagnetic field. Such a K-drive will in general 
give rise to a non-zero bulk acceleration field A = "VyV on Mi. In principl^ 
one could maintain any prescribed state of acceleration Aq with V = Wq by the 
addition of a further drive j™'''^* of some nature, provided 

> / A iT^\ rin,elastr^ i , rin,EAl ^ . rin.extr^ i 

Ck huAo{K) = fj^' [J:t]+fK [St], 

where j™'*^'"^** isolates possible drives in the medium due to internal 
non-electromagnetically induced stresses and friction with the environment. In 
particular if the medium is maintained in a state of uniform rotation in the ab- 
sence of interaction with electromagnetic fields by mechanical torques then the 

''This would be difficult in practise if electromagnetic stresses in a deformablc medium led 
to significant strains producing a dynamic change in its shape or volume. 
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presence of any electromagnetic induced torques can (in principle) be compen- 
sated by additional mechanical torques in order to maintain a state of uniform 
rotation. Using the orthogonal decompositions with respect to U above and the 

1 , . - 7 in^EM 7 jin,EM . r in. EM i • , x j • j- i- 

reiation ifjciTjf — dJj^ + LjjPk , such an mtegrated mstantaneous 
external X-drive may be written: 

U,in,extrr^ i _ I ■ ^ in, EM _ I r U,in,EM I jU,in,EM 



Ik [^t\ = ~ / iudTj^ = CuPk ^ Jr > (26) 



where Stokes' theorem has been used in the last term. The precise nature of 

in"' ' 



this compensating drive will depend of course not only on the form of j-™'^*^ 



but also on the electromagnetic constitutive relation for the medium. 

For unbounded media the notion of a total integrated drive may not be 
meaningful. In these circumstances one may be able to deduce contributions to 
isT-drive pressures from the last term in or by integrating discontinuities 
of over surfaces instead of integrating contributions to ijjdTj^ over volumes. 
This follows simply from the bounded medium case above by writing 

jU,in.EM I jU,out,EM , / \ i tU\ 

Jr - -Jr + ^12{Jr), 



where the discontinuity A12 is across the boundary of the history in spacetime. 
Both integrals now require a knowledge of the pull-back to the interface of J^°^* 
in the vacuum region. However one may relax the condition that the medium be 
bounded and assert that the instantaneous _fC-drive exerted on any finite spatial 
volume Vt in an infinite history having an interface 5*4 = dVt with the rest of 
the medium and the vacuum, by drive 3-forms that contain discontinuities over 
part of St, is given by 

nU,in,ext\^f 1 I f U.in,EAl I jU.out^EM I \ f jU \ /or^A 

Jr [Vt\ = - / I-uPr - Jr - / ^12(Jr)- (^7) 

JVt -JSt -JSt 

In static situations the first integral on the right is zero and in some cases the 
second integral can be evaluated in terms of fields or their sources in the vacuum 
region. 



4 The Abraham and Symmetrized Minkowski 
KiUing-Drives in Non- Accelerating Media 

As noted in the introduction there have been a number of different proposals for 
the electromagnetic stress-energy-momentum tensor. In the tables in Appendix 
A, the properties of a number of these [27] : are displayed using the notation 
introduced above. For a recent derivation of these tensors from a variational 
approach see [281 . They permit a ready evaluation of the iiT-drives on media 
with bulk 4-velocity field V and specified constitutive properties. For conve- 
nience, separate entries for the Abraham tensor have been provided for media 
at rest relative to a frame U (i.e. V — U). To calculate forces and torques 
it is only necessary to recognize that they are defined by particular spacelike 
Killing vectors K. So in an inertial frame U with U{K) = many entries will 
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simplify. Separate entries are also given for situations when K is timelike and 
equal to U since this facilitates computation of energies and powers. In the 
calculation of the divergence of the Abraham stress-energy-momentum we have 
made explicit those terms that depend on the bulk acceleration field A = Vy V 
of the medium. Such terms are of course absent for media at rest or moving 
with arbitrary linear velocity in an inertial frame (VJ7 = 0). It is only in spe- 
cial circumstances that there are simple relations between the Minkowski and 
Abraham tensors. These relations arise mainly for simple media at rest in an 
inertial frame. Thus if 

= eoC-e'^, b*^ = AioA'rh^ (28) 

in a stationary homogeneous medium (with constant scalars and /i^), then for 
all K, the are identical for rSM^^A/ ^^^^ j,ab (ggg ^^^Ae in Appendix A). 
However the corresponding are all different. In problems with time- harmonic 
electromagnetic fields the Cup^ have zero averages over time, implying that the 
corresponding time-averaged if-drives ( ) are the same. 

In a simple non-accelerating {A — 0) medium that may be inhomogeneous 
there is a simple relation between and t^*^ (see table in Appendix A) for 
arbitrary time-dependent fields. Given a spacelike translational Killing vector 
K one sees from the table that (since U (K) ~ 0) 

= r|*^-i.^(if).f/, 

where the 1-form 

= * (^^e^ A A C/ - cod^ A A f/j . 

Using the simple constitutive relations (PSI) this 1-form can be expressed in 
terms of the spatial fields e'^, h'^ and the form U 

= -(l-AA2)*fe^Ah^AC/V (29) 
Co ^ ^ 

where Af = y/f-rfJ-r is the refractive index of the medium. Thus 

s''{K)irU^~-(l-N^)iu fe^Ah^AC/Aif) = - (l - N^) fe^Ah^Aif' 
Co V / Co V 

and so 

= + ^{M'-l){e'' Ah^ Ak). (30) 

By contracting the exterior derivative of this expression with U — -^dt (defining 
a non-accelerating reference frame) one deduces a relation between the 3-form 
force density j^]^'^ = iudr^ associated with the Abraham drive 3-form, and 
■pU,SM _ ^^^^s^M associated with the symmetrized Minkowski drive 3-form: 

^u,A ^ ^^^^^M + _L£,^((AA2_i)e^^h^)Ai?, (31) 
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using 



d(t> = d(t>~UhCu(t> = d(t> + dt A Cd^4>, (32) 

and noting that (W^ - l) {e^ Ah" AK^ is a spatial S-fornfl. The second term 

on the right in ((3T|) has been termed the Abraham force. In static field situations 
it is zero and the two force densities are equal. This is also the situation in the 
presence of time-periodic fields if one takes the time-average of this equation 
over a time period. For smooth pulsed fields the second term will in general 
yield a finite contribution to a total force when is integrated over a finite 
time interval. 

Returning to the constitutive case for an arbitrarily moving medium () one 
can use the results in the previous section to facilitate a comparison between the 
Newtonian force density and that based on the choice with j-^'^*^ = ^k^^ ■ 
Thus one must calculate iu drf/^^ for a Killing field K that generates spatial 
translations and express the result in terms of the polarization and magneti- 
zation 1-forms p'^ and respectively. A local constitutive relation between 
G and F will induce a local constitutive relation between and m'^ and the 
electromagnetic fields and b*^ and the permittivity and permeability of the 
medium. Since dF — 0, 

2dT^^' = -{diKGA*F-iKGAdi.F-FAdiK*G). 

Writing G — e^F + 11, one has 

diKG A-kF ^ CkG a *F - iKdG A-kF ^ F A -kCxG - iKd{eoF + A-kF 
= F ACk*G - ixdU A*F = F A ixd *G + F A dix *G- F A i^ixdU 
= F Aixj + F Adix *G - F A ★i^dH. 

Hence 

2dr|^'^ = ~F AiKj + F A-kiKdlV + ixG Ad-kF. 

Since 

d-kG = eod * F + d-kll = j, or d-kF= — {j — dirU), 



and 



ixG Ad-kF = -G Aixd-kF = -—GAix {j ~d-kU) 

eo 

= -— G A ixj + —G A iKd-kU = -F A ixj - — A ixj + G A ixd-k — , 
eo eg eo eo 



one has 



2dT|*^ = - {2F +—] Aixj + F A ^^ixdll + G A i^d * — . 

eo / eo 



^If one replaces U hy V with V^V^ ^ in II29II one is dealing with accelerating media and 
the Abraham drive J-"^' then acquires additional terms in JSTJ. Such additional terms are 
often overlooked when considering the effects of the so called Abraham force (or torque) . This 
may lead to inconsistencies in processes involving accelerating media. 
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This can be decomposed in terms of e*^, b^, p^, m^, and real and induced 
polarization and magnetization sources yielding^ 

-2d4'' = [pl^e^iK) + i^J^ A ^ + i(p^ + p'Dv^iK)) A U 

+ [poi-u{J" + Jp ) A K + 2p^e^(A:) + A k + ^i^T> J^m A «: 

-copZh'^iK)] AU- 2coB^ A ixp"" - tV^'" A iKp"" 
-iKP^Jn A + coJA'B^ A + ^iif A p^, 
where M*^ = #Af mC/. Contracting with C/ yields 

jr^'^*^ = -^c;dr|*^ =P^e^(7^) + t^J^ A n + ip,^e^(A-) 

+ i^g7^p'' A «; + ai^i^T, A k - icop^b^(A-) (33) 
+ 2^(P'' + Pp + ^Moi^ C^'' + Jp^) A K. 

The first two terms on the right constitute the Lorentz force density due to any 
real charge fP and electric current in the medium. The third and fourth 
terms constitute the Lorentz force density due to induced polarization charge 
Pp and electric current . These are precisely the force densities that arise 
in the non-relativistic electrically neutral smoothed model ([CT) above based on 
the motion of point charges. The fifth and sixth terms constitute the Lorentz 
force density due to magnetization charge and magnetization current J^. 
These do not arise in the Newtonian model above since the point charges 
in that model were not endowed with intrinsic magnetic moments. 



5 Conclusions 

The general theory of drive-forms has been developed starting from the van- 
ishing divergence of a total stress-energy-momentum in an arbitrary spacetime 
using the language of exterior systems. It has been shown that a decomposition 
of a drive 3-form, relative to a unit timelike vector field, on a spacetime with suf- 
ficient Killing vectors yields forms that can be associated with force and torque 
densities in continuous media. For material subject to the electromagnetic con- 
stitutive relations proposed by Minkowski it has also been shown how the com- 
putation of the electromagnetic force density on a non-accelerating polarizable 
perfectly insulating medium in an electromagnetic field can be effectively carried 
out in terms of a particular split of the total stress-energy-momentum into parts 
describing its inertial and electromagnetic properties. In paper II this theory is 
applied to homogeneous and inhomogeneous dielectric media where it is argued 
that a means of discriminating between a split into the electromagnetic stress- 
energy-momentum tensors proposed by Abraham and Minkowski (and possibly 
others) can be explored experimentally using rotating media. 
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A Summary of Electromagnetic Stress-Energy- Momentum Tensors 





MINKOWSKI 


SYMMETRIZED MINKOWSKI 


rpEM 


iaF ®i'^G+\-k{F A *G)g 


T^M = K (g) i°G + iaG (g> i°F + *(F A *G)g) 




{iRF A*G - F AiK * G) 


-if- (ikG A*F- F AiK*G) 


V • 

rpEM 


i d * (F A *G) - ( Vjf „ G) {Ff) - G( V^F) 


i ( d* (F A *G) - (VxaGK^F) - {^x^F){FG) - G(V^F) - F(V^)) 


Jk 


-Cif \e^(K)DU + h^(K)BU - ■^U(K)e^ A 
L ^0 

_i(et'.d^+b^-h^')«] 


- (e^- d^ +hU.hU)K- U{K) (^e^ A h^' + cod" A b^)] 


P^K 


Ck [cod^ Ah^ AK-^ 

+ '^U{K) {e^'-d'-' + W'- W') #lj 


CiL [(-LeC^ Aht^+cod^^ Ab'^) Ai^-L + [/(ii:) (et^.dt^+bt^-hC^)#l] 




-^e^' A li" 

CO 


1 (-Let^Ah^'+cod^'Ab^) 






KeC^-dC^ + b^^-hC^)*! 





ABRAHAM {U 7^ V) 


ABRAHAM {U = V) 


rpEM 


1 ( io-F <8> i"G + iaG (g) i°F + *(F A *G)s - ® F - y ® s^) 


5 ( iaJ' ® i"G + iaG ® i°F + ★(F A *G)g - s'^ (g) - C/ (g) s^') 




-Sf (iKGA*F-FAix*G-s^(/f)*y- y(if)*s^ ) 


-if (ikG A ★F - F A ix * G - s'^(if ) * t/ - U(K) * s*^ ) 


V ■ 

rpEM 


V . T^M - 1 ((V • s^)y + (V • V)8^ ) - ^iV) - F{G(A)) + G{fCa)) 


V . T'S'^ - 1 ({V • s'^)U + (V ■ U)a'^ ) - *([/) 




-if. [eU{K)T)^ + 6P{K)-EP + h^(i<:)B^ + b^(i<')H^ 

- {e^' ■ d*^' + b*^' • h*^') K - U{K) (^e'^ A h'^ + cod^ A b*^) 
-iu [k Mv* s^) + 2V{K)iu * s^] 


-Cif [| (e^(A:)Dt^ + d^(K)E'^ + h^{K)BU + h^{K)H^) 

— 2\ "tD ■n — — U )G AO J 




if- [(^e^ Ah'^ + cod'^ Ab^) AK-^ -U(V)U(K)*s^ 

+U(K) (e*^ • d^^ + b^^ • h'^) #1 + U(y)iu * A K + V{K)Uu ★ s^] 


Ck [^et^ A A i^-L + iC/(is:) (e^.d^ + b*^ ■ h*^) #l] 




i (^e*^ A h*^' + cod*^' A b^^) + ^Uuiv * + C/(V)i[/ * 


CO 


pYj 


i (e^^ • d^ + b^^ • hf^) #1 - U{V)Uu * 8^ 


i(e^'.d^+b^'.ht')#l 



In these tables we have introduced the 1-form: 

= i,{ivF Aiv*G AV ~ivG Aiv*F AV) 

where the spatial fields e^,b^,d^,h^ are defined by the orthogonal splits of 
F and G relative to the 4-velocity field V of the medium. Furthermore, we 
have defined H'^ = #h'^. The vector fields U and V on spacetime are timelike, 
unit, normalized with a metric g of signature (— 1,+1,+1,+1) so that U{U) = 
V{V) = —1, and describe the state of motion of the observer and medium 
respectively. It is also convenient to introduce, in terms of the vector field U 
and any vector field Y on spacetime the projection = Y + U{Y)U, which 
induces the projection operator on p- forms Hjj — l + U Aijj- The spatial Hodge 
map # with respect to U is defined by ★! — U A^fl, and maps spatial p- forms to 
spatial (3— p)-forms. For any Killing vector field K on spacetime (i.e. Cxg = 0) 
we denote the 2-form by k. If iiT is a Killing vector then so is K where 
^0 is any constant. Hence the physical dimensions of quantities that depend on 
the Killing vector field K will depend on the physical dimensions of f o • 

In the table we have isolated those terms in the divergence that depend 
explicitly on the bulk medium acceleration A = VyV. Among the remaining 
terms there is a set that can be expressed in terms of a (2,0) -tensor field 5* on 
spacetime defined, for any vector fields X,Y, by 'if{X,Y) = (^(X))(F) where 

^{X) = {\/xF){i^) + F{i^^G)-iVxG){i7F)-G{i^^F)-^{£xF){i^) 
-lF{{£xg-^){ixG)) - \F(i7c^G) + i(/:xG)(jxF) + \G{{Cxg-^)(ixF)) 

+ \G{i7c^F) + l{Cxg){F{GG)) - \{Cxg){G{i^F)) 

with F{Y) = iyF and G{Y) = iyG. If {Xa] denotes an orthonormal frame 
with dual cobasis {e''} one has g^^ — ry"^ Xa <E) Xi,, g = ryah (E) e** where 77"'' 
and rjab are both diagonal matrices with 77'^'^ = 7700 = — 1 and rjij = i]^^ = Sij 
for i,j = 1, 2, 3. The contraction operator ix^ is abbreviated ia and i'^ = ■q"''' ib- 
The spacetime Hodge map is denoted * and the canonical 4-form on spacetime 
is ★! = e" A A A in an orthonormal basis and is -^/jdet g\ d'^x in any 
coordinate system {x"}. 

B Notation 

The natural mathematical language to discuss the differential properties of ten- 
sor fields on spacetime and their relation to integrals over material domains is 
in terms of differential forms and their associated exterior calculus [18] . In this 
section a brief summary is given of the relevant notation used in subsequent 
sections. A key concept throughout involves the role of the spacetime metric 
tensor field and possible isometrics that it may possess. A spacetime metric 
tensor field 5 is a symmetric bilinear form on spacetime that can always be 
represented in a local cobasis of differential 1-forms {e"^} as 

g = -gO (g) (g) (g) (8) e^. (34) 
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If {Xt,} is the dual local basis of vector fields on spacetime defined so that 
e°'{Xi)) = (a, b, = 0, 1, 2, 3) one has the induced contravariant metric 

5"^ = -Xo0Xo + Xi0Xi+X2 0X2+^3 0X3. (35) 

If a is any given l-form, one has an associated vector field a defined so that 
/3(a) = g~^{a,l3) for all l-forms (3. Since is symmetric this will be abbre- 
viated a = g~^{a). In a similar way if X is any given vector field one has an 
associated l-form X = g{X). 

In the £(-orthonormal basis {e"} one has a canonical local 4-form denoted *1 
and defined to be e° A A A e^. The Hodge map ★ induced by ★I maps p- forms 
to (4— p)-forms on spacetime (Benn & Tucker 1988; see also Appendix A of paper 
I). The metric also uniquely defines the covariant derivative Vx with respect to 
any vector field X. This has the property Vxg = and VxY — VyX = [X, Y] 
for all vector fields X, Y. In this expression [X, Y] denotes the commutator 
bracket. While Vx has a type-preserving action on any tensor field the exterior 
derivative d is defined to act only on antisymmetric tensor fields (differential 
forms) and has the property d o d = 0. Contraction of any p-form /3 with X is 
denoted ix(3- The interior operator ix is a graded derivation defined so that 

ixiaAp) = {ixa)Ap+{-lfaAixP, (36) 

for any p-form a and g-form /3. If p = 1, one defines ixct = a{X), and if p = 0, 
ixct = 0, for all vector fields X. One has the useful relations 

★ = (_i)P+i$ 
ix*^ = -k{^AX) 
Cx^ = ixd^ + dix^, 

for any j>form $ on spacetime where jCx denotes Lie differentiation with respect 
to X. 

The spacetime divergence operator V- takes a simple form if one uses the 
(7-orthonormal basis above. Acting on a symmetric covariant tensor T it defines 
by contraction on the first argument the l-form 

3 

V-T = ^(Vx„T)(X«,-) (37) 

a=0 

where {X°-} = {— Xq, Xi, X2, X3}. For a symmetric tensor T it is sufficient 
to write the right hand side as (Vjc^ T){X'^). A local spacetime isometry with 
respect to is a local diffeomorphism that preserves this metric. A vector field 
K that generates such a diffeomorphism is called a Killing vector field and it 
satisfies Lxg = in terms of the operation of Lie differentiation with respect 
to K. The operators ■k,d,ix,'Vx,^x offer a powerful computational tool-kit 
when working with differential forms. 

For any smooth p-form $ in a bounded regular region of a manifold one 
can express the integral of over A4 in terms of the integral of $ over the 
boundary dM of M: 

[ = f (38) 

Jm JdM 
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This is a statement of Stokes' theorem for p-forms. 

The Gibbs calculus of vector fields in 3-dimensional Euclidean space is read- 
ily exposed by correspondences induced by the exterior operations above. A 
space of 3-dimcnsions may be considered as a particular hypersurface in space- 
time. If the metric above induces the metric 

g = ig) + (g) + (g) (39) 

on this hypersurface it is Euclidean and one may introduce the Euclidean canon- 
ical form #1 = A A e'^ by restriction. Since spacetime is assumed time- 
oriented one may employ a futijxe-pointing timelike unit vector field U on space- 
time (with g{U, U) = —1 and U = e^) to fix a coherent orientation by relating 
#1 to 7^1 by 

★ 1 = UA#1. (40) 

If the hypersurface is given as t = constant for some time; coordinate t, an 
inertial frame exists in Minkowski spacetime with U = -^-^ = —Xq such that 
V [/ = 0. Throughout this article the constant cg denotes the speed of light in 
the vacuum. Any p-form a on spacetime is termed spatial with respect to such 
a U if iua = 0. An over-dot will denote (Lie) differentiation with respect to the 
coordinate t so 

a = cqCuQ- = jCg^a. (41) 

On a Euclidean hypersurface in spacetime, exterior differentiation of spatial 
p-forms (f> is denoted d(f> such that 

# = - U A jCu'P = d(j) + dt A jCdt4>- (42) 

One has the following relations 

#1 = -i^i,i = -i,u 
= cj> 

for all spatial p-forms ^. Furthermore if v, w denote Euclidean vector fields 
in the Gibbs notation corresponding to the vector fields v, w for some spatial 
1-forms V, w then 

div V corresponds to #d#v ; (43) 

curl V corresponds to #dv ; (44) 

V X w corresponds to #(w A w) ; (45) 

grad tp corresponds to dtjj, (46) 

for any scalar ^ field on spacetime. 
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C Electromagnetic Fields on Spacetime 



We suppose that a polarizable material continuum is given in terms of a set of 
piecewise smooth material properties that determine its interaction with classi- 
cal gravitational and electromagnetic fields. The classical macroscopic Maxwell 
system for the electromagnetic 2-form F on spacetime can be written as (Tucker 
& Walton 2009): 

dF = and d*G=j, (47) 

where the excitation 2-form G depends on the interaction with the medium and 
the 3-form electric 4-current j encodes the electric charge and current source^. 
Such an electric 4-current describes both (mobile) electric charge and effective 
(Ohmic) currents in a conducting medium. To close this system in a background 
gravitational field, electromagnetic constitutive relations relating G and j to F 
are necessary. 

The history of a particular observer field in spacetime is associated with an 
arbitrary unit future-pointing timelike 4-velocity vector field U. The field U 
may be used to describe an observer frame on spacetime and its integral curves 
model idealized observers. An orthogonal decomposition of F with respect to 
any observer field U gives rise to a pair of spatial 1-forms on spacetime. The 
1-form spatial electric field and 1-form spatial magnetic induction field 
associated with F are defined with respect to an observer field U by 

= iuF and cob^ ^ iu * F. (48) 

Since g{U, U) = -I and iue^ = iuh^ = 0: 

F = AU -*{coh^ AU). (49) 

Likewise the 1-form spatial displacement field and the 1-form spatial mag- 
netic field associated with G are defined with respect to U by 

d'-' ^ iuG and — = ic/ ★ G, (50) 
Co 



so 

G = d" AU AU j , (51) 

with ijjd^ = iijii^ = 0. At the history of any sharp interface between different 
media, given as the piecewise smooth (non-null) spacetime hypersurface / = 0, 
the system of Maxwell equations is supplemented by interface conditions on the 
fields F and G 



[F] 



Adf = 

(52) 



All electromagnetic tensors in this article have dimensions constructed from the SI di- 
mensions [M], [L], [T], [Q] where [Q] has the unit of the Coulomb in this system. We adopt 
[g] = [L'^] , [G] = [j] = [Q] , [F] = where the permittivity of free space eo has the dimen- 
sions [Q^T^ M^^ L~'^] and cq = ^g^^^g denotes the speed of light in vacuo. Note that the 

operators d and V preserve the physical dimensions of tensor fields but with [g] = [L^], for 
p-forms a in 4 dimensions, one has [*a] = [«][L^~-^f]. 
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where [H] denotes the discontinuity in the field H across the hypersurface. The 
3-forni js on the hypersurface is non zero if it supports a real current 3-form 
there. 



D Time-Dependent Maxwell Systems in 3-Space 

The spacetime description above is natural for the Maxwell system since it makes 
no reference to any particular frame in spacetime. However to make contact with 
descriptions in particular frames or non-relativistic formulations a reduction in 
terms of frame dependent fields becomes mandatory. The spacetime Maxwell 
system can now be reduced to a family of parameterized exterior systems on M.^. 
Each member is an exterior system involving forms on depending paramet- 
rically on some time coordinate t associated with U. Let the (3 + 1) split of the 
electric 4-current 3-form with respect to a foliation of spacetime by spacelike 
hypersurfaces with constant t be 

j = — A[/ + p^, (53) 
Co 

with iijj^ — ijjp^ = and = p^^l, where ^p^ are the spatial electric 
current density 2-form and spatial electric charge density 3-form respectively. 
The differential operator d on spacetime forms is adapted to those spacetimes 
(such as Minkowski spacetime where gravity is absent) that can be foliated by 
hypersurfaces with constant coordinate t where U — j^dt- Then, from (|Tfl) . 

dj = (54) 

yields 

dj^ + p^ = 0. (55) 

The (3-1-1) split of the spacetime covariant Maxwell equations (|T7)) with respect 
to U = —c^dt becomes 

de^ = -B^ (56) 
d'B" = (57) 

dW^ = J^' + t>^ (58) 
dD^ = p^ (59) 

where = #d'^ and B*^ = #b'^. All p-forms {p > 0) in these equations are 
independent of dt, but have components that may depend parametrically on t. 



E Electromagnetic Constitutive Relations 

The two 2-forms F and G in the macroscopic Maxwell equations on spacetime 
are fundamentally related by smoothing the microscopic sources of the elec- 
tromagnetic fields in the medium. In many circumstances one then relies on 
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phenomenological relations for closure relations. In such relations the excita- 
tion form G is in general a functional (possibly non-local in spacetime) of the 
Maxwell form _F, its covariant differentials, thermodynamic properties, defor- 
mation properties and the state of motion of the medium: 

G = Z[F,VF...]. 

Such a functional may induce non-linear and non-local relations between d'^ , h'^ 
and e'^, b^. Elcctrostriction and magnetostriction arise from the dependence of 
Z on the deformation tensor of the medium and its covariant derivatives. For 
general linear continua, a knowledge of a collection of constitutive tensor fields 
Z^"^") on spacetime may suffice so that 

N 

G = ^^('''[V'^i^,...]. 

r=0 

In idealized (non-dispersive) simple continua, one adopts the idealized local re- 
lation 



G = Z{F), 

for some degree 4 constitutive tensor field Z, parameterized by scalars that 
depend on the medium. In the vacuum G — e^F where eo is the constant 
permittivity of the vacuum. Regular lossless, non-conducting, linear isotropic 
media can be described by a bulk 4-velocity field V of the medium, a real 
relative permittivity scalar field > and a real relative permeability scalar 
field fir > 0. In this case, the structure of the tensor Z follows from 

G = eoerivFAV-—*(iv*FAv) 

IXr ^ J 

1 \ ^ (60) 

Cr I ivF hV + —F. 

fir J fir 

For inhomogeneous media the relative permittivity and permeability scalars tr 
and fir will not be constants. In a general frame U comoving with the medium 
{U = V),^ yields 

d^ — eofire^ and = (^oMr)~^b^, (61) 

which are the familiar closure relations for simple (idealized) electrically neutral 
isotropic non-dispersive polarizable media. 




E.l Minkowski Constitutive Relations for Moving Media 

The above algebraic constitutive relation (pO)) involves the bulk 4-velocity field 
of a simple medium. It is straightforward to find the induced relations between 
the fields {e'-^, b'^, d*^, h^} relative to an observer in a frame U . From 
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ivF = {V)U - U{V)e'^ - * cpb^ ^U ^V 



iyi^F^ ★(e^ AU AV)+ cob^(y){7 - CQU{V)h^ 



ivG = {V)U - U{V)d'^ - ★ 



AU hV 



Co 



h^(^)~ C/(V^) 



and ([?T|) . one has 

cob^ = 
d^ - 

Co 

Inserting these in the constitutive relations (j6ip yields 

d^(y)f7 - U{V)6" - ★ A (7 A i/) = e [ e'^'(l^)(7 - (7(y)e'^ - * (coh^ A A V^)] 

Co * (d^ A t/ A i/) + h"{V)U - U(V)h" = /i-i [ * A f/ A i/) + b^(l/)[/ - U(V)h"]^ , 

where e = e,eo,/^ = /^rMo in terms of the relative permittivity scalar and 
relative permeability /i^. Contracting with U gives the relations 



U - 



Co 



Co 



h^(y)=At-ib^(T/), 



(62) 



which yields 



U{V)d^ 



AU AV\ = e 



Co 



f/(F)b^ - Tk-I^AC/AFj U{y)W^ - Co ★ (d^ A [/ A 



(63) 



A laboratory Minkowski inertial frame is described by C/ = — 9t in inertial 
Cartesian coordinates {t,x^,x'^,x^} with {e° = —codt,e^ = dx^,e^ = dx^,e^ = 
dx^} and, relative to U, the medium 4- velocity V has the orthogonal decompo- 
sition: 

Co 



V ^ j[U 



where the spatial field is the Newtonian velocity fielj"! of the medium 

relative to U, 7"^ ee J(l - 4) and z/^ = g(V^, V^). With these definitions it 
V '^o 

follows that: 

g{U,V) = -j and U AV = —U AW^. 

Co 

Using these in (j63p and rearranging yields 



7^ ( [/ A A ^ I = e 



b*^ + ★ A A 



★ ( t/ A A b'^ 



[/ A A d^ 



(64) 



^^In inertial Cartesian coordinates = V]"^ , f-^ (t, x^, x'^, x^) t;^. 
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or using the spatial Hodge map # defined by U: 



+ # A 



- # (V^ A 



(65) 



These are the constitutive relations first written by Minkowski. For some pur- 
poses it is useful to decouple these expressions and express e'^ and h*^ directly 
in terms of d*^, b'^ and V*^. 



E.2 Minkowski Constitutive Relations for e^(d^, b^, V^), h^(d^, b^, V^) 

Taking the exterior product of ([55)1 with V'^ yields, 

/xV^ A h*^ = A b^ - V'^ A # (y^ A + fxV^ A # (v^ A d*^ 

For any spatial 1-form with respect to U one has: 

A # (V^ A a'^) = -V^ A ★ ([/ A A a'^) = - * i^u (u AY'^ A a*^) 
= i.2^(c7Aa^) -a^(V^)*(c7AV^) 

Thus 

tiY^Ah^ = V^Ab^ + 4#e^-^^^#V^-M^^'#d^ + Md^(V^)#V^ 



(66) 



= A b^ + ^#e^ - Aij.2#d^ _ 



using the identity d'^(V'^) = ee^(V^), obtained by contracting (US]) with V*^. 
Since ## = 1 



/i# ( V'^ A h'^ 



# ( A b^) + ^e^ - fii^'d^ - -^{1 ~ eMC^)d^(V^)V^. 



Substituting this into the first relation of (j65p yields 



e^CQ 2 ' 



Mc^ 1-^ d^-(.MC^-l) #(v^'Ab^) 



d^(V^), 



Since 



eg 



J\f- 



'2 ' 



where Af^ — erfir is the square of the medium refractive index scalar, this may 
be written 



AA2-r^)e^ = :^(l-r^)d^-(AA2-l) (#{V^Ab^)- 



d^(V^) 
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Furthermore from (|65p 

= # (V^ A d^) - + ei^'b^ - ^ {eficl - l) b^(V^)V^, 

V / Cq /iCq 

using b^(V'^) = /ih^(V'^), obtained by contracting dHS]) with V^. Substituting 
this into the second relation of (|55|) yields 

Thus, the constitutive relations (l60l) can also be written 



^ - 9 f 9 (-# r -■') - 

In the non-relativistic limit (to first order in — ) these constitutive relations 
become 



eoer V <^ 



(67) 



Mo Mr \ fi r Mr 

F Polarization and Magnetization 

The polarization 2-form 11 in spacetime is defined by 

n = G-eoF. (68) 
The second macroscopic Maxwell equation may then be written 
eod-kF = j - d*U = j + jp, 

where 

jp = -d-kU (69) 

will be called the electric polarization current 3-form. With respect to any 
observer frame U its orthogonal decomposition is 

n = p^Ac7-*f— Ac7^ =p^'ac/--m^, (70) 

V Co / Co 

where M'^ = #m'^ and we call 

= ij/II and =i(7*n 

Co 
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the spatial polarization 1-form and magnetization 1-forni respectively relative 
to U. The Hodge dual of n has the decomposition 

*n = *(p^AC/) + A [/ = + -m^ A [/, (71) 

Co Co 

where = . From dCT]), dH]), dMl) and ([701) it follows 

d^ = eoe^+p^ and h^^o^b^ + m^. (72) 

From dSni), dZOl) and ^ one finds 

j„ = -dpC^-^A[7 = -dP^ + C/A£c;P^-^A[7 
Co Co 

= -dP^ + — (cqCuP'^ - dm^) A U 

Co 



Co V / 

Writing the orthogonal decomposition of jp with respect to U as 

Co 

it follows that 



iujp^^(p'^-dm'^^ and = -(i^/ ★ jp) * [/ = -dP^ 



Co ■ Co 



In the frame U , Jp and denote the induced electric polarization current 
density spatial 2-form and induced polarization charge density spatial 3-form 
respectively. In a similar manner 

dn = dp'^AC/-— dM^^dp'^Af/-— dM^ + — t/A/Ic/M^ 

Co Co Co 

= — (codp'^ + /:c/M^) A (7 - — dM^. 

Co Co 

with the orthogonal decomposition j„i = dll = A [/ + where 

— = -^uhn = - f codp^ + -M ^) and = - * J™) * C/ = - -dM^ 
Co Co V Co / Co 

denote the induced magnetization charge current density spatial 2-form and in- 
duced magnetization charge density spatial 3-form respectively in terms of p*^ 
and M*^. 



G Derivation of the Model Force Density 3-form 

The key steps in deriving the expression for the force density 3-form in §2 of 
paper I are outlined (equation 2.14), with the aid of the identities 

#aA/? = (73) 
CpOt = ij^da + dij^a (Cartan's identity) (74) 
d(aA#/3) = (75) 
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for any spatial l-forms a,/?. From equation (2.8) of §2, paper I, we have 
Tk = i[d(e^-p^ + b^-m^)+#/:a.(p^Ab^)] A^, 



(76) 



with K = =ffK. The first term may be expressed 

d(e^-p^) = <M~,e^ ^ C~,e^ - i~,de^ 

using (j73l) and (l74t . The last term can be written 

(p^Ab^) = #(p^Ab^+p^Ab^) =#(p^Ab^-p^A#de^ 

= #(p^Ab^) + #(#de^Ap^), 

using (j56l) . Finally, the second term is 

d(b^.m^) = dz^m^ = C^m^ - i^dm^ 

= /:^m^-##z^rfm^=/:^m^-#(#dm^Ab^), 

using (j73l) and (17^ . From these expressions 

r^e^ + r^m^ + # ((p^ - #dm^) A b^)" 



2J"if 



A K 



£^e^ + £^,m^ + # (#Jp^ A b^^ 



A K 



A K, 



using (|73l) and (1751) . Using Cartan's identity in the first term here gives 



— di^r A k) — e'^ A di^K — A i^dn 
= dZpT, (e^(if)#l)-e^'Ad#(i^Ap^) 
= d(e^(X)P^) +e^Ad#(p^ Ai^), 

using ([74| and ((75|) and the fact that K — and so dK = 0. In the last term 
one has 



(77) 



A d#(p^ A ^) = A dixV^ = A Ck^^ - A i^dP^ 
= e^A#/:Kp^-e^(i^)dP^ 
= £KP^AE^ + p^e^(i^), 

since X is Killing and using ([75|. This yields 

p;;'e^(/0 + /:/fP^ A + d (e^(AOP^) 

and so the force drive 3-form becomes 

= ]-\p^B"{K) + igjJ^ Ah + CkT?^ A¥F +Cgjxa^ AK + d{e^{K)I' 



C-r,e^ A K 



29 



H Orthogonal Decompositions used in Deriving 
the Drive 3-form Associated with the Sym- 
metrized Minkowski Electromagnetic Stress- 
Energy-Momentum Tensor 

In this section, the key steps used in deriving equation (4.6) of §4, paper I are 
outUned. For U{K) = 0, one calculates 

— AiKj = UAi—iKJ^AU + iKP^^ 



eo V^o 



1 fpt/ A [/ - -mA a (-ikJ"" AU + iKp"") 
eo V Co / V^o / 

— AUA iKp" - A ixJ" AU - —M" A ixp" 

eo eoCfl eoCo 

( — A iKp^ - MoM^ A IkjA AU- — A iKp^ 
\eo / eoCo 

-p^{K)p^ + fioi-vJ^ Ak) AU- —M^ A ikP^ 



,eo / eoCo 

FAiKj = (e^ AU-cqB'^^ A(^^iKj" AU + iKp'^^ 

= e'^ AU A iKP^ - A ixj'^ AU- cqB'^ A ixp^ 
= (e^ A iKp'^ -B^ A iKJ^) AU - cqB'^ A ikP^ 

= (e^{K)p" +igiJ^ Anj AU-CoB^ AiKp^ 

F A MKdU = ixdn A*F= (^^ixJ^ AU + ixpZ^ A (e^ + cqW^ A U^ 

= —IkJ^ AU AYF + iKp^ A + CQiKp^ Ah^ AU 
Co 

-iKJii A E^ + cqIkpZ, a ) a t/ + ikpZ, a 

cq / 



--i^J^ AK + Coplh^{K)\ AU + iKpl A 

Co ^ / 



GAixd-k— = ——iKGAd-k'n. = —iKFAd-k'n. — iKllAd-k — 

eo eo eo 

-ixFAd^U = -(^iKe^ AU-coixB^) a(^-^J^ AU-p^ 

= iKe^ AUAp^- IkB^ AJ^ AU- cqZxB^ A 

= - (e^(/f)p^ + ixB^ A J^) AU- CoiifB^ A 

= (-e^ (i^ - A k) A [/ - Coi^B^ A 
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eo eo V Co y V cq ^ J 

= -iKp" AU Ap^ - -^IkM" AJ^ AV - — i^M^ A 
eo ^ eoc^ eoCo ^ 

= - i-^^{K)p'i + MoixM^ A jf\ AU- — z^M^ A 

Veo ^ J eoCo ^ 

= ('-lpt'(i^)p^ - txoi^J^ Ak)aU- — ZkM^ a p^. 

These expressions may be combined to produce the required 4-form as 
described in §4, paper I. 
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